Skyrmions and edge spin excitations in quantum Hall droplets 
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We present an analysis of spin-textures in Quantum Hall 
droplets, for filling factors v ~ 1. Analytical wavefunctions 
with well defined quantum numbers are given for the low-lying 
states of the system which result to be either bulk skyrmions 
or edge spin excitations. We compute dispersion relations 
and study how skyrmions become ground states of the Quan- 
tum Hall droplet at v ^ 1. A Hartree-Fock approximation is 
recovered and discussed for those spin textures. 

PACS number: 73.40.Hm 

Spin excitations of a two-dimensional (2D) electron 
gas in the filling factor v — \ regime are the subject of 
increasing interest due to the experimental observation 
[0-^ of skyrmions. They have been analyzed both 

by a classical nonlinear-cr model (CNLSM) ||,|^ and by 
a Hartree-Fock (HF) approximation. These descrip- 
tions ignore quantum fluctuations; furthermore quantum 
numbers as modulus (5^) and third component {Sz) of 
the spin, and total [M) and center-of-mass [Mcm) an- 
gular momenta are not well described. A superposition 
of mean-field wavefunctions has been proposed [|lO[ which 
have well defined M and but still neither Mcm nor 
are represented adequately. Current descriptions of spin- 
textures have always been considered in infinite [||"0J^] 
or periodic (sphere) p|,[7|jlC|] systems. Therefore, they can 
not be directly used to analyze available experimental in- 
formation for small systems. |TT| - [l3| 

In this letter we address the two previously mentioned 
restrictions of present theories: 

1) We analyze spin-textures in a system of N electrons 
moving in 2D, in the presence of both a high perpendicu- 
lar magnetic field B and a confinement potential (Quan- 
tum HaU droplets (QHD)). 

2) We obtain analytical many-body wavefunctions of 
spin excitations with all quantum numbers (M, Mcm, 
and Sz) properly described, i.e. we present a microscopic 
description including quantum fluctuations. 

We also obtain the energies of these excitations, which 
turn out to be low- lying excitations &t v = 1, some of 
them becoming the ground state (GS) close to = 1. 

We consider interacting electrons moving in the xy 
plane, confined by a parabolic potential (characterised by 
a bare frequency wq). We assume that B is high enough 
as to project the Hamiltonian H onto the lowest Landau 
level with both spin up and down. In symmetric gauge. 
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tion operators, respectively, rrii are single-particle (SP) 
angular momenta and aj label spins. We consider ma- 
trix elements Vm-im^m^nn for the Coulomb interaction, 
although we have checked that our conclusions do not 
depend on the details of the electron-electron repulsion. 
The single particle part H^^ has a spectrum 
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where lUc = eB/{m.*c), fl = [uj^ + 4uJq]^^^, g is the Lande 
factor and /is the Bohr magneton. Due to the symme- 
try of the system the spectrum separates in subspaces 
with well defined M, Mcm, S^,Sz- The energy ordering 
of states within a subspace does not depend neither on 
luq nor on B; only the relative positioning of energies in 
different subspaces do. 

Some of these subspaces have dimension 1, so the state 
which spans each one of them is an eigenstate of the 
many-body Hamiltonian. This is the case for single Slater 
determinants comprised of the k lowest m SP states with 
spin I and the N — k lowest m SP states with spin |, 
where k = 0, l,...iV/2. We label these eigenstates as 
|CAr-fe)i and we refer to them as "compact" states, due 
to the the fact that they have minimum M compatible 
with a given 5*^ = {N-2k)/2. Numerical calculations Q 
show that compact states are GS for some ranges of the 
parameters defining the QHD. For example, on increasing 
B the transition from the paramagnetic IC*^^!) — 2) to 
the ferromagnetic \C^) (i' = 1) occurs by flipping spins 
through all |C^_j,). However, these are not the only 
GS in the process, other GS with small Sz appearing 
in between any spin-flip transition. Small-^z states 
appear as wefl as GS for < 1. 

We want to describe these highly correlated low-S'^ GS 
and, more generaly, low-lying excitations close to = 1. 
It has been shown that there is a branch of low-lying 
excitations at v = 1 which resemble magnetoexcitons of 
|C'5v)- 10 analogy, in order to describe spin-charge 
excitations, we try the ansatz of magnon-like excita- 
tions: OljC%^^), with =EmOmcJ„±A,iC,„J. It 

is known that low-lying excitations have Mcm = 0. |lj] 
As for charge excitations dX v = 1, |jl^ imposing this 
condition fixes the values of the coefficients Om- We find 
|]l9| several possibilities for 0\^^ (with A > 0): 
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where and c are the electron creation and annihila- 
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It can be shown that (Oj.^)", n being an integer and 

Oj.^ any of the previous operators, also generate states 
with well defined Sz, M, and Mcm- For that purpose, 
these operators can act on any state |C^_j.), except 

(A|^)", which requires A < fc (note that therefore none 
of the (A^)" act on the = 1 state). However, (Oj_^)" 
generate states without a well defined S"^ . Due to the 
Zeeman term, the lowest energy state in a (M, 5'^) sub- 
space has maximum possible Sz- Let us call V the oper- 
ator which projects onto the subspace — Sz{Sz + !)■ 
po| So we generate, from compact states, sets of states 

^(^±A)"I^Af-fc) which have all quantum numbers prop- 
erly defined. In fact, the use of the V operator opens 
another possibility for O^-type operators: 
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acting on any state |C^_j,) docs not have well defined 

Mcm, but P(sl)"|C^_j,) has Mcm = for all n. 

In order to check whether these ansatz states are eigen- 
states of Hamiltonian (|l|), we have performed numeri- 
cal diagonalisations for < 9. The overlap between 
7^(C^^^)"|C^_^) and the exact lowest energy state in 
the {M,Sz) subspace they belong to is always higher 
than 0.99. This is so for O being T,A, or S, and 
for all A, 71, and k which satisfy the restrictions stated 
above on the applicability to these operators. There 
is only one case in which two different states of the 
form ViOl^)"\C%_f^) belong to the same subspace: 

r{A\)"\C^_^) and ■P(SI)"|C^_;,) Wc find that, for 
A: > 0, the former one has always a high overlap with the 
exact lowest energy state of the subspace it belongs to, 
while the later one (when orthogonalised to the former 
one) has an overlap again larger than 0.99 with one of 
the low-lying excited states in that subspace. For fc = 
(i.e. — 1), there are not states generated by A|, and 
7'(s|)"|C^_j,) becomes the lowest energy state of the 
subspace with M = A^(A^ -l)/2 + n and Sz = N/2 - n. 

Now we investigate the physical nature of the different 
kinds of states. Let us focus first on states of the form 
7^(A|)"|C^_j^). Simpler wavefunctions can be obtained 
by the inverse transformation to the one which projects 
a skyrmionic BCS-like single Slater determinant into a 
state with well defined M. We obtain 
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being the vacuum state. From this form, |5'((/3, ^)) 
is readily recognised as a HF wavefunction describing 
skyrmions with topological charge 1; it depends on 

two parameters: is a variable fixing a broken-symmetry 
direction of the spin in the xy plane, while parameter ^ 
controls the skyrmion size (i.e. its M and Sz)- From 
eq. (0), S'f^(r) - (r2 - 2^^)/{r^ + 2^^) at large dis- 
tances, as for skyrmions in the CNLSM, which implies 
a size Q A = ^/V2. Remarkably, such dependence of 
Sf^{r) holds even in the case of skyrmion sizes compa- 
rable with the magnetic length Ib — \/h/m*Q,^ regime 
well beyond the range of applicability of the classical 
model. A technical point, of relevance in the interpre- 
tation of these states, is that in eq. (|[) the projector 
V is not included. Otherwise \^{if,^)) would not be a 
BCS-like single Slater determinant. Corrections intro- 
duced by V are proportional to S*^^, so for n << N and 
A^ >> 1 they are negligible. In the small A^ limit, or for 
n ^ N/2, the skyrmion is "deformed" and its identifica- 
tion as such becomes fuzzier. For instance, we find (for 

N even) V{AI)W')-^\C],_,)=V{tI,)'^/^\C%). 

By a similar argument, we identify excitations created 
by A^ as skyrmions with topological charge A. Recall 
that, although skyrmions with A > 1 are never the lowest 
energy excitation branch in this model, they are eigen- 
states appearing from partially-ferromagnetic compact 
systems with, at least, A electrons with spin |. 

A physical identification of states generated by e| or 
is not so direct. From the form of the operators 
both of them are edge spin-textures. Unlike what occurs 
for skyrmions where 5*^ changes sign as a function of 
r, their Sz{r) is always positive. They resemble spin- 
waves in that they appear as excitations of any compact 
(partially- ferromagnetic) system. 

In an infinite system, electron-hole symmetry can be 
invoked to obtain other spin-textures. In particular 
antiskyrmions arise as (Aj^)-type excitations of states 
|C^^°°) of holes. For the antiskyrmion with topolog- 
ical charge —1 and Sz = n, we obtain the wavefunction 
(E™=i"^"^/^4^_i,|C™,|)"cn,T|C^^^). For a finite sys- 
tem, the hole state complementary to the electron state 
co,i\C%_^_i) is not a compact state and it has Mcm 0. 
In this case, antiskyrmions cannot be built up by opera- 
tors of the form Oj_^. 



Numerical diagonalisations for small A^ show that, on 

to 



increasing B, the spin-flip transition from |C^^^ 
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occurs through a ladder of states of the form 
7'(Al)"|C^+\,^,P and V{TUr\C^^_^). In this letter 

we concentrate on the transition from |C^_i) to \C%), 
which is sketched in Fig. 1. Energies of skyrmion- like 
states = -p(A|)"|CjV_i) take the form 

(^^SK\^SK) {M,Sz) + -to + aAn+fJA\^2 J 

where Eq is the interaction energy of the initial compact 
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^) in this case. Both and /3a have ana- 



lytical expressions in terms of V, 
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which can be 

calculated, by using Wick's theorem. We find Q!a < 
and /3a > 0, so skyrmions in a QHD behave as confined 
bosons interacting via a two-body repulsive interaction. 
Figure 2 shows the dispersion relation of the skyrmion 
branch for N = 30. Energies for 7'(tLi)"|C^) have the 
form of eq. (^) , with different expressions for , ax and 
f3r- We find that /3t is negligible for all N. The transi- 
tion from \Cl^_-^) to |C^) presents two cases: 

i) for N < Nq, where A'o depends on the Zeeman en- 
ergy but is Nq < 10, the transition is through a lad- 
der of states with decreasing first, until it reaches its 
minimun possible value and increasing Sz afterwards, as 
depicted in Fig. 1. 

ii) for N > No the minimum of the skyrmion branch 
is not at the end of the branch (n = N/2 — 1), as shown 
in Fig. 2 for = 30. The transition is through all the 
skyrmion states with n less than a certain n^^ , which 
depends on g, going directly to |C^) afterwards without 
passing through any of the 'P(tL-^)"|C5^) states, n^^ 
labels the spin of the largest skyrmion which becomes 
GS when varying B. 

The spin of the GS is a measurable quantity which may 
be controlled varying the Zeeman energy by means of a 
tilted B. Figure 3 shows n^^ as a function of Zeeman 
energy, calculated from ^'f ^) and the HF approximation 
(HF energies can be obtained using eq. (g)). The HF 
approach gives slightly smaller values for n^^ . This is 
so because, although the HF dispersion relation is very 
accurate (see inset of Fig. 3(a) and note the small energy 
scale) its minimum is shifted to smaller n. 

Experimental information on the transition from 
\Cjf_i) to \C%) is already available: single electron ca- 
pacitance pl| and transport | 13| e xperiments in quan- 
tum dots have been interpreted | |l3| , p^ in terms of highly 
correlated GS's appearing in numerical diagonalizations 
performed for N < 5. Here, we have identified those 
states as skyrmions and edge-spin excitations and shown 
that skyrmions appear for all N. 

To complete our discussion on skyrmions we analyse 
the thermodinamic limit. In the g = limit, CNLSM 
predicts that skyrmions in an infinite system have smaller 
interaction energy than \Clf^^) by an amount AEint — 

-y^T:/32{e'^/4nelB)- [|| Figure 4 shows n^^ at g = and 
the interaction energy of the skyrmion of size n^^ as a 
function of N. n^^ grows slowly with N, while AEint 
tends to the classical value. 

From numerical diagonalisations we know that, for 
small N, ||l5|,|l9| only two kind of excitations become 
the new GS of the system when increasing B: charge 
magnetoexcitons JI^\C%) ||1^ or edge spin excitations 
7^(s|)"|C5^). However, for an infinite 2D system, an- 
tiskyrmions must be the new GS of the system for v ^ I. 
Since we do not have analytical expresions for wavefunc- 
tions and energies of antiskyrmions in the intermediate 
N regime, we cannot conclude whether they are the new 



GS in this case. With this caveat in mind, we concentrate 
on charge and edge-spin excitations at u = 1. For arbi- 
trary N, expressions for magnetoexciton energies have 
been reported. 0211 ^e find that edge spin excitation 
energies are of the form (||) with negligible P^. There- 
fore, their energies are linear in n, i.e. in M, see Fig. 
2 (notice that this is equivalent to a dispersion relation 
which is quadratic in linear momentum /c, as occurs for 
excitations produced by a broken continuous symmetry). 
Which type of excitations has lower energy depends on 
both N and Zeeman energy. In the g = limit, we 
find that the edge spin branch has lower energy than the 
charge excitation branch for N < Nmm and vice versa 
for N > N,nin where Nmin — 100 for the Coulomb repul- 
sion. On increasing Zeeman energy, the charge branch is 
favoured. For instance, for GaAs {g = 0.44), Nmin ~ 20 
for tiLUQ — 2meV while for hujQ = bAmeV the charge 
branch has lower energy than the spin branch for any 
number of electrons {Nmin = 0). 
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FIG. 1. Sketch of M and of the different states which, 
in some conditions, become GS of the QHD around v = 1. 
Compact (•), skyrmions (o), edge spin excitations (x), (o) 
and charge excitations (A) are shown. 
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FIG. 2. Dispersion relations of skyrmions (Ai), edge spin 
excitations (T_i), (E) and charge excitations (J) for N = 30. 

AM is the excess angular momentum with respect \Clr_i). 
Parameters for E^^ {M,Sz) are taken such that \C%){*) and 
have the same total energy. 



FIG. 3. Exact (continous line) and HF (dotted line) Sz for 
the largest skyrmion which becomes GS when varying B (in 
the text refered to as n'^^) as a function of g for iV = 30 and 
(a) TiijJo = 2meV, {h)hu}o = 5.4meV. Inset in panel (a) shows 
the comparison between exact (dots) and HF (dotted line) to- 
tal energy for skyrmions (in units of /(AnelB))- Parameters 
for E^^{M,Sz) are taken as in figure 2. 



FIG. 4. n (dotted line) and interaction energy of the 
skyrmion of Sz = n^^ with respect to that of |C^_i) (conti- 
nous line) as a function of N, in the g — limit. * marks the 
interaction energy for skyrmions given by the CNLSM in an 
infinite system. All energies are in units of /{AirelB)- 
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